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FOREWORD 


The  present  work  arose  out  of  a  need  to  consider  sound  absorption  processes 
in  viscoelastic  aster ia Is  other  than  naturally  occurring  dissipation  by 
aolecular  relaxation*  The  latter  becomes  especially  weak  at  low  frequencies  and 
should  be  supplemented  by  alternative  mechanisms.  Of  several  possible 
electnwechanical  and  magnetomechanical  dissipation  processes  we  choose  to 
investigate  absorption  due  to  mechanically  induced  eddy  currents. 

The  work  was  carried  out  in  the  Nonmetallic  Materials  Branch  of  the 
Materials  Division  of  the  Research  and  Technology  Department  under  an  1R  grant. 
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CHAPTER  1 

ACOUSTIC  REFLECTION  FROM  A  UNIFORMLY  MAGNETIZED  VISCOELASTIC  PLATE 
ON  A  CONDUCTING  GROUND  PLANE 


The  purpose  of  this  report  is  to  assess  the  performance  of  an  acoustic 
attenuator  in  the  form  of  a  viscoelastic  sheet  attached  to  an  electrically 
conducting  backing,  the  sheet  consisting  of  either  alternating  magnetized  and 
unmagnetized  strips  of  material  or  a  distribution  of  magnetized  disks  in  an 
unmagnetized  lattice.  In  the  present  paper  we  consider  the  situation  where  the 
magnetized  regions,  assusmd  to  be  electrically  nonconducting,  are  well 
separated,  thus  focussing  attention  on  absorption  due  to  induced  eddy  currents 
in  the  backing  or  ground  plane  while  neglecting  interactions  between  magnetized 
regions.  We  also  treat  the  magnetized  regions  as  if  they  were  infinite  in 
extent  as  far  as  the  acoustic  response  is  concerned,  thus  neglecting  edge 
effects  due  to  differences  in  the  elastic  responses  of  magnetized  and 
unmagnetized  regions.  These  elastic  responses  may  be  different  not  only  because 
of  the  magnetic  moment  but  also  because  of  basic  differences  in  the  elastic 
moduli  and  densities  of  the  materials. 

The  two  simplest  magnetic  structures  in  this  category  are  uniformly 
magnetized,  either  in  the  plane  of  the  sheet  or  perpendicular  to  it.  In  the 
present  paper  we  focus  on  uniform,  perpendicular  magnetization. 

In  order  to  come  up  with  a  readily  calculable  problem,  we  consider  only 
one-diamnsional  geometries.  We  realize  at  the  outset  however,  that  finite 
uniformly  magnetized  plates  do  not  generate  uniform  magnetic  fields.  We 
therefore  have  to  make  a  number  of  simplifying  assumptions  to  reduce  the 
magnetoelastic  interaction  problem  to  a  one-dimensional  one.  Consequently,  the 
resulting  theoretical  developments  will  be  highly  approximate.  In  addition,  we 
have  to  make  certain  mathematical  approximations.  Nevertheless  it  is  hoped  that 
the  results  will  convey  at  least  the  order  of  magnitude  of  the  expected 
attenuation. 

Consider  a  uniformly  magnetized  plate  with  the  magnetization  vector 
perpendicular  to  a  conducting  ground  plane  and  bonded  laterally  to  an 
unmagnetized  viscoelastic  material,  as  shown  in  Figure  1.  Let  the  shape  of  the 
plate  be  either  a  circular  disk  or  a  long,  rectilinear  slab.  We  will  argue  that 
we  may  replace  the  magnetization  by  an  equivalent  surface  current  distribution, 
as  shown  in  Figure  1.  We  may  then  solve  for  time  varying  fields  above  and  below 
the  ground  plane  as  functions  of  time  varying  surface  currents,  I Q(z,t)  (ampere- 
turn/  cm),  where  the  time  variation  of  10  is  ultimately  due  to  acoustic 
vibrations  in  the  plate. 
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Ht  itiuae  that  the  local  magnetic  dipole  moment  is  frozen  into  the 
viscoelastic  material.  Hence* it  does  not  change  because  of  the  relative 
displacement  of  two  neighboring  points  of  the  material.  The  only  change  in  the 
magnetisation  envisioned  here  is  that  due  to  changes  in  the  local  density  of 
djipple  moments  which  occurs  as  a  consequence  of  such  local  relative  displace¬ 
ments  of  the  material.  The  magnetization,  Mg,  is  the  magnetic  moment  per  unit 
volume: 


“  dV 

In  line  with  our  goal  of  considering  a  one  dimensional  problem,  we  assume 
that  the  motion  of  the  plate  in  response  to  a  normally  incident  acoustic  wave  is 
strictly  in  the  s-direction  as  shown  in  Figure  1.  This  makes  Mg  and  10 
functions  of  z  only.  Furthermore,  in  this  case  7  x  M0  *  0  and  since  the 
SKMit  general  source  function  for  the  vector  potential,  aQ,  is: 


Vo 

4» 


F  X  Mg 


y*r 


we  find  that  it  is  a  function  of  the  surface  current  I„  ■  H„  xi  only  -  even 
when  a  compressional  wave  exists  in  the  slab.  Throughout  the  motion  of  the  slab 
we  may  therefore  calculate  the  vector  potential  on  the  basis  of  the  current 
source  distribution  shown  in  Figure  1.  The  source  potential  becomes: 


a/ 


where  I  (z0,t)  is  a  function  of  the  local  displacement  u(zQ,t).  Thus 
coupling  between  the  magnetic  and  the  viscoelastic  problem  occurs  because  of  the 
presence  of  u  (zQ,t)  in  a0. 

The  dependence  of  I  on  u  may  be  calculated  as  follows.  Suppose  we  are 
given  a  distribution  f(s)  on  a  <  z  <  b.  Consider  the  new  distribution  ht(z')  on 
«'t  <  *'  £b't»  «*t  ^  b't  ^b,  which  is  obtained  from  the  coordinate  trans¬ 

formation  s'c  ■  g(z,  t),  where  we  indicate  that  the  transformation  may  be  a 
function  of  a  parameter  t,  such  as  time.  Since: 


h(z')  dz' 


b 

■  y  hlg(z.t)]  U  dz 


we  have 
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h[g(z,t)]  ■  f  (z) 


(5) 


In  the  present  problem  we  are  working  with  the  transformation: 


z'(t)  ■  z  +  u(z,t)  z  g(z,t)  (6) 

and  the  distribution  f(z)  •  IQ  ■  constant  on  the  interval  0  £  z  £  L,  which 
exists  when  no  acoustic  disturbance  is  present  in  the  plate.  Identifying 
h(z'(t))  with  l(z,t)  at  any  given  time,  we  find: 


Hence,  in  the  presence  of  a  compressional  wave,  the  current  at  z'  *  z  +  u(z,t)  is 
modified  by  9u(z,t)/az  at  z.  For  the  long  wavelength  limit,  which  is  the  one  we 
are  interested  in  here  3u/>z  is  small  and  we  therefore  disregard  this  dependence 
and  replace  I(z')  by  IQ. 

In  order  to  solve  for  the  acoustic  response,  we  need  the  magnetic  volume 
force  exerted  by  the  total  magnetic  field  at  any  point  inside  the  slab.  Ue  are 
going  to  assume  that  this  force  may  be  calculated  in  a  certain  average  sense 
from  the  same  equivalent  surface  current,  I0,  that  is  being  used  to  calculate 
the  magnetic  field.  Suppose  we  have  solved  the  magnetodynamic  boundary  value 
problem  for  the  current  I0(z,t)  and  have  found  the  vector  potential  at  z  £  0, 
A(z,t).  Let  A(z,w)  be  the  Fourier  transform.  The  resulting  Lorentz  force  on 
an  equivalent  current  loop  will  generally  not  be  normal  to  the  plane  of  the 
slab.  But  in  keeping  with  our  aim  of  simplifying  this  problem  to  the  point 
where  only  vertical  forces  and  motions  remain,  we  consider  only  the  y-component 
of  the  total  magnetic  field.  We  also  consider  only  the  case  of  a  rectilinear 
strip  of  material  of  width  26  and  thickness  L,  since  this  geometry  is  of 
primary  interest  to  us.  The  case  of  the  circular  disk  can  be  analyzed  in 
exactly  the  same  way.  We  give  the  corresponding  vector  potential  below.  Hence 
the  Fourier  transform  of  the  Lorentz  force  per  unit  height  and  unit  length  of 
slab  on  e.g«  a  right  hand  surface  current  element  of  strength  I(z,w)  at  z  is: 


F(z,w)  ■  I(z,u)  X  By(z,w) 


-  -I(z,w) 


»A(z,m) 

iZ 


(8) 


/.v.v, 
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In  order  to  obtain  an  expression  for  the  total  magnetic  volume  force  (force/unit 
volume)  normal  to  the  plane  of  this  current,  we  first  add  the  force  on  the  left 
hand  current  filament.  This  force  is  in  the  same  direction,  even  though  the 
field  is  in  the  opposite  direction,  because  the  current  is  reversed.  Dividing 
by  the  width  of  the  strip,  2A,  we  obtain  the  average  magnetic  volume  force  or, 
store  accurately,  we  obtain  a  simple  analytical  quantity  which  we  will  use  in 
place  of  the  true  average  force  inside  the  material: 


Pm(z,»>  -  -  ~  [l(z,«)  X  “”7*— J  (9) 

From  this  force  we  have  to  substract  the  equilibrium  magnetic  force  which  exists 
when  no  acoustic  field  is  present  because  the  displacement,  u(z,t)  in  the 
solution  of  the  acoustic  problem  represents  the  displacement  from  position 
coordinates  reached  not  only  under  the  influence  of  static,  mechanical  forces 
but  also  under  the  influence  of  magnetic  forces  due  to  the  static  magnetization, 
Mq.  Hence  the  resultant  volume  force  is: 


Pn(z,u) 


Io  r»A(z,tn)  >Aq(z)1 
a  [  az  ■  a*  J 


UO) 


This  is  the  Fourier  transform  of  the  local  magnetic  force  due  to  all  sources; 
internal,  due  to  other  magnetic  moments,  as  well  as  external,  due  to  induced 
eddy  currents  in  the  ground  plane.  The  latter  force  will  enter  the  dynamic 
equations  explicitly,  whereas  the  magnetic  force  due  to  other  internal  moments 
will  be  taken  care  of  by  means  of  a  phenomenological  modification  of  Hooke's 
Law.  However,  for  the  moment  we  will  not  distinguish  between  these  various 
sources  and  calculate  the  total  vector  potential  resulting  from  all  time  varying 
internal  (equivalent  surface)  and  external  (eddy)  currents.  The  part  of  the 
resultant  potential  that  is  due  to  internal  sources  will  be  discarded  later. 

Let  the  applied  vector  potential  for  unit  stationary  surface  current  at 
z  ■  zQ,  be  a0(z,z0).  This  yields  the  slab's  static  magnetization  which  is 
unmodified  by  the  presence  of  the  ground  plane  if  we  assume  that  the  magnetic 
permeability  of  the  ground  plane  is  that  of  free  space,  which  we  do.  Then  the 
applied  potential  due  to  a  suface  current  element  to  ID  at  z  *  zQ  ♦  u(zQ,t)  at 
time  t  is  a0(z,zQ  +  u(z0,t)).  Expanding  this  to  first  order  with  respect  to  zQ, 
or  more  accurately,  taking  the  Frechet  differential  with  respect  to  u(zQ,t)  at 
z  »  zQ,  we  write: 


9a0 

a^(z,Zg  +  u(Zg,t))  m  Bq( Z,Zq)  ♦  u(zg,t)  j 

°  u“o 


■  a0(z,z0)  ♦  Aa0(z,z0) 


(11) 
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Mow  let  a  and  Aa  be  Che  Fourier  transforms  of  Che  solutions  for  the  total 
vector  potentials  corresponding  to  aQ  and  bag,  respectively.  Then: 


AqU) 


a(z,*0)  dze 


,.>  -  / 


[a(z,zQ)  +  Aa(z,z0,u>)J  dzQ 


Hence,  the  total  local  magnetic  pressure  is: 


Pm(z,u)  -  -  y*  Aa(z,z0,u»)  dze 


aa(z,z0,m) 


a  r  .  .  >«vz,z, 

—  /  jj- 

O 


where  we  recall  that  the  vector  potential  in  Equation  (14)  is  understood  to  be 
only  that  part  of  the  total  potential  Chat  is  due  to  induced  eddy  currents  in 
the  ground  plane. 

This  pressure  has  to  be  accounted  for  in  the  dynamic  equations.  When 
Pm  •  0  these  are: 

Hooke's  Law: 

r  i2  >2u  aP  fl 

>Icl  teK’-it  11 


Newton's  2nd  Law: 


a^u  _  ap 


where  P  is  the  mechanical  pressure,  p  the  mass  density  and  [c]  the  speed  of 
sound  in  the  absence  of  magnetic  forces.  In  order  to  obtain  the  corresponding 
tnagnetoviscoeiastodynamic  equations  we  simply  add  the  magnetic  volume  force  on 


•  7  .*  /  / 
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the  righthand  tide  of  Equation  (15b)  and  modify  Hooke's  Law  implicitly  by 
changing  the  wave  speed  in  the  material  from  [c]  to  c*c(Mq).  In  other  words, 
the  internal  magnetic  interactions  change  the  elastic  properties  of  the  material, 
conceivably  even  its  viscoelastic  properties,  by  moving  the  equilibrium  state  of 
the  dynamic  response  to  a  different  point  on  a  nonlinear  response  characteristic. 
To  account  for  this  change,  we  simply  modify  Hooke's  Law  in  the  way  indicated 
above.  The  resulting  equations  are: 


2  _  a* 

pc  hit  at 


(16a) 


u  _  _  ap 
,t2  "  '  az  + 


(16b) 


Ue  regard  c  as  an  emperically  determined  quantity.  It's  dependence  on  MQ 
should  be  kept  in  mind  throughout  the  remainder  of  this  paper.  When  the  change 
i,n  performance  of  the  magnetic  slab  is  investigated  as  a  function  of  the 
interaction  with  the  ground  plane,  via  Pm,  it  must  be  kept  in  mind  that  as 
PB  (specifically  I0)  changes,  so  does  c(M0),  unless  we  assume  that  the 
longitudinal  modulus  of  the  material  in  which  the  magnetic  moments  are  embedded 
is  adjusted  to  conpensate  for  the  effects  of  internal  magnetic  interactions. 
Conversely,  if  we  change  c  while  keeping  Pa  fixed,  we  will  ultimately 
encounter  incompatible  combinations  of  the  two  and  may  obtain  unphysical  results. 

Since  u(z,t)  appears  under  the  integral  sign  on  the  right  hand  side  of 
Equation  (16b)  this  is  an  integro-partial  differential  equation.  It  can  be 
solved  readily  by  means  of  perturbation  techniques  and  we  indicate  below  how 
this  is  done.  A  quicker  way  is  by  resorting  to  the  following  approximation.  We 
replace  u(z0,t)  in  Pm  by  its  average  value  ono<  i0<l: 


u(z0,w) 


1 

L 


L  _ 

J  u(z0,w)  dz0  s  u(z0,u)  s  a(w) 
o 


(17) 


and  obtain: 
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In  this  fora  the  equations  nay  be  integrated,  readily.  We  are  still  taking 
account  of  the  coupling  between  the  acoustic  and  electromagnetic  fields,  but  are 
sacrificing  the  point  by  point  amplitude  and  phase  relationship  that  exists 
between  the  local  displacement  and  the  y-component  of  the  magnetic  field.  In 
view  of  other  approximations  that  were  made  earlier,  this  seems  to  be  a 
reasonable  step  at  this  point.  The  solution  for  u(z0,t)  will  be  a  function  of 
e  and  later  on  we  will  have  to  solve  for  a  using  Equation  (17).  Since 
a(«)  will  occur  linearly  on  both  sides  of  this  equation,  this  step  is 
trivial. 

For  the  reflection  coefficient  we  need  the  steady  state  solution  of 
Equation  (16)  with  PB  as  defined  by  Equation  (18).  Like  <*(«),  the  magnetic 
fields  in  Equftion  (18)  remain  unspecified  for  the  moment.  Introducing  the  time 
dependence  e""1***-  and  setting  v  *  Ju/*t  we  obtain: 


pc 


2  4L 

dz 


i«P 


(19a) 


-iwpv 


(19b) 


Substituting  Equation  (19a)  into  Equation  (19b)  one  finds: 


d*v 

dz2 


♦ 


k2v 


(20) 


This  inhomogeneous ,  second  order,  ordinary  differential  equation  must  now  be 
solved  subject  to  certain  boundary  conditions.  Introducing  the  incident  and 
reflected  waves: 


*o 

vo 


AlP  -  -  BlP‘"o 


-ika(z-L)  +  B  _ik„(z-L) 
z-L) 


-i- 

Po«o  1 

L 


-  BlPlk° 


(z-L)j| 


and  requiring  that  both  v  and  P  be  continuous  across  the  boundaries,  we  find: 


z-L  ;  Poccv(L)  •  Bj  -  Ai  -  P„c0v0 


(21) 


z-L 


iPc  dv 
k  dz 


Ai  +  Bi 


*o 


(22) 


i  v  jTT V V  v  V  •  i  *  • 


*•''•,***  V*" »  •  V* 

?  •  •  .*  •  *  »  -•  ’  •  •  • 


and: 


m?  I1  P  q  U  ww»r  row  »■■  I'vr 
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Defining  Che  known  input  admittance  of  Che  backing  as  G(0)*-v(0)/P(0) ,  one  has: 


*  -  0  ;  G<0)  ~  -  v(0)  -  0 


(23) 


Combining  Equations  (21)  and  (22)  we  obtain  the  boundary  condition  on  v  at  z*L: 

2ikAi  dv  /  ip0c0k 

z  *  L  ;  -  -  - - I - 

pc  dz  \  pc 

and  from  Equation  (23)  at  z*0  we  have: 


(24) 


0 


dv  /  ik  \ 

"  dz  +  \pc  G(0)  )  v 


(25) 


Suppose  the  solution  of  Equation  (20)  has  been  found.  The  final  step  is  to 
calculate  the  reflection  coefficient,  *  I Bx/Aj_J  ^ .  Using  Equations  (21)  and 
(22)  one  finds  that: 


p0c0v(L) 

P(L) 


V-a  .  _ 

ikfc£o\ 

v(L) 

l  ♦  a 

\  PC  / 

dvt 

dz|E«L 

»  m 

Poco  G(L) 


(26) 


where,  by  definition,  G(L)  is  the  input  admittance  of  the  magnetic  layer. 
Equation  (26)  will  therefore  yield  the  answer  to  the  problem  once  v(z)  has  been 
found. 


The  most  compact  expression  for  v(z)  is  obtained  with  the  aid  of  the 
Green's  function,  g,  of  this  problem.  This  function  oust  satisfy  mixed 
homogeneous  boundary  conditions  according  to  Equations  (24)  and  (25),  i.e.,  for 

g  we  must  have: 


z  -  0 


Di(g)  i  ?g 


♦ 


ik 

pc  G(0) 


8 


0 


(27) 


z  *  L 


02(g)  *  ?g 


ikp0c0 

pc 


0 


(28) 


Referring  to  the  general  solution  of  this  type  of  problem  in  three  dimensions  as 
given  in  standard  texts,  the  solution  of  the  present  problem  may  be  written: 
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i,«)  -  [g(z,0.w)  -  v(z,«)  -•-•^■° -~j  |  ^  _  q 

ik  > 

-  —  j  Pm(z0,w)  g(z,z0,w)  dz0  ;  0  <  z  <  L 


where  g  is  the  solution  of: 


dg2(z.s0iw)  >» 

- -  ♦  kg(z,z0,w)  ■  -  «(z-zc)  (30) 

In  Equation  (29)  the  usual  surface  integral  has  been  reduced  to  function 
evaluations  at  the  two  boundary  points  and  the  volume  integral  has  become  a  line 
integral.  Using  the  boundary  conditions,  Equations  (24),  (25),  (27),  and  (28), 
one  finds: 


v(z,m)  •  +  *2 


where: 


2ikAx 


g(z,L,«) 


-*/ 


pm<zo«“)  g(z.z0><»>  dzQ 


The  function  therefore  yields  the  reflection  coefficient  of  the 
nonmagnetic  layers,  whereas  ^  yields  the  effects  due  to  magnetization. 

Had  we  not  used  the  average  displacement,  a(w)  in  Equation  (18),  we 
would  now  have  a  double  integral  on  the  right  hand  side  of  Equation  (33)  with  v 
appearing  in  the  integrand.  Equation  (31)  could, there fore, be  used  as  the 
starting  equation  for  a  perturbation  expansion  of  the  solution.  For  instance, 
with  the  aid  of  Equations  (14),  (31),  (32),  and  (33)  one  finds  that  to  first 
order,  the  solution  for  v  is: 
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v(t,«) 


2ikAt 

pc 


g(z,L,l) 


g(z0,L,w) 


* 


3a(z' ,z0,w) 
3z '  3ze 


g(z,z' ,«)  dz' 


(34) 


with  remainder: 


Av(z,w) 


/ 


/ 


L  L 

f  f  v(tQ,w) 
o  o 


*a(t' ,t0,w) 
3t'  3t0 


*  gU0. *'»*»> 


3a(z 1 ,z0,w) 
3z '  3z0 


g(z,z',u)  dz' 


dz0  dt'  dt0 


(35) 


and  8  follows  from  Equation  (26).  Since  a  (z,z0>u)  is  proportional  to  u0I0 
and  g(z,z0,w)  to  1/k,  it  is  clear  that  the  expansion  parameter  in  this 
perturbation  series  is: 


*  •  (fe) 

The  ratio  of  the  magnetic  field  energy  density,  u0I02  and  the 
mechanical  kinetic  energy  density  pc2  in  Equation  (36)  is  a  small  quantity 
in  all  real  magnetoelastic  materials.  Furthermore,  it  may  be  assumed  that 
1 4k I  is  of  order  unity  or  larger  at  all  frequencies.  Hence,  X  is  a  small 
number  and  the  expansion  converges  rapidly.  On  the  other  hand,  the  smallness  of 
this  parameter  also  tells  us  that  the  effects  we  are  going  to  observe  will  be 
weak  since  the  solution  with  finite  Pm  will  differ  little  from  that  for  a 
nonmagnetic  slab. 

In  order  to  find  v  explicitly,  we  need  g.  By  definition: 


g 


-1 

*<>(yi.y2) 


yi(z)  y2U0J  1  2  <  zo 

y2(*)  yi(*oM  *  >  *o 


(37) 


where  y^,  and  y2  are  two  linerarly  independent  solutions  of: 
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and  40,(yiy2)  is  the  Wronskian: 

Ao(yi.y2>  •  (yi  Jp  -  /2 

Evidently,  g  will  satisfy  its  boundary  conditions  if  yi  and  y2  are  chosen 
such  that: 

z  ■  0  ;  Dl(yi)  “  0 

s  ■  L  ;  D2(y2)  “  0 

We  use  the  two  linearly  independent  functions: 

y^  ■  sin(kz  -  e^) 

and 


y2  *  ainlk(z-L)  -  «2] 
with  Wronskian: 

Ao<yi.y2>  “  "k  ”  *1  ”  *2) 

Applying  the  boundary  conditions  yields: 
tan[*i)  *  -ipc  G(0) 

and: 


tan[«2]  •  - 


IOC 


Aoco 


for  the  two  integration  constants  in  Equation  (31)  and  the  problem  is  solved 


12 
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For  g(z,zto,w)  one  finds: 


1  Tsinllcz-ei]  sin[k(L-z0)  - 

-  1  L  k  -  *1  -  »2>  J 

Z  + - ►  ZQ 


Z  <.  Zo 


(46) 


Z  >  Z0 


Removing  the  singularity  In  the  slope  of  the  Green's  function  from  the  Integral 
in  Equation  (33)  leads  to  a  more  explicit  form  of  the  solution: 


z 

v(z)  -  C1  yi(z)  -  ^  f 


[yi (z)  y2(z0)  ~  y2<z>  yi(zo>l  Pm^o* 


dzr 


where: 


‘1  ‘  c2  *  K  / 
o 


Pm(*0)  y2<zo)  dzo 


and: 


I 

2ikAi  ( y2(L)\ 

C*’-  PC  \  t  ) 

Observing  that  the  Wronskian  is  independent  of  z,  we  note  that: 


^ yi .y2)  |L "  (y*y2’  “  ^i’*  |  L  *  -^i'  ’ 1  yi)  |  L 


-y2(L)  D2[yi(L)] 


Substituting  into  C2,  we  find  that: 


2ikAx 

C2  nr  n  fv_fl.ll 


(47) 


(48) 


(49) 


(50) 
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Before  writing  down  the  reflection  coefficient,  we  calculate  a(w). 
Writing  Equation  (31)  in  the  form: 


v(z)  ■  Aif}(z)  ♦  a (»)  f2U) 


(51) 


where: 


fl(*)  *  g(*,L,«) 


(52) 


end: 


g(z,z0,w)  dz0 


end  using  Equation  (17)  one  finds: 


(53) 


®(w) 


l-»2 


(54) 


where: 


and: 


(55) 


(56) 


For  the  purpose  of  estimating  the  performance  difference  between  magnetized 
and  unaagnetized  states  we  write  B  in  the  form  8  *  80  ♦  48,  where  80  is  the 
reflection  coefficient  in  the  absence  of  magnetization  and  48  is  the  part 
depending  on  Mg.  A  measure  of  the  effectiveness  of  the  magnetization  for  two 
otherwise  similar  materials  is  given  by: 


14 
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R  *  [ »a « 2  -  I a0 1 2 ]  ■  2  Real  [48  80*]  +  I4BI2 


We  therefore  want  to  calculate  both  8  and  80.  With  the  aid  of  Equations  (24), 
(26),  (51)  and  (54)  we  find: 


1-8  _  _/t v 

i  ♦  a  **oco  g(l)  *  -  p0c0 


Setting  ♦j  *  0  *  ?,  w  obtain: 


(  f i(L)  (1  -  *2)  +  4l  ^2^^  ) 

| [2  +  Pqco  fj(L)l  (1  “  42>  +  p°c°  *1  f2(L)  ) 


1-8 


Vo  G(L)P 


=  pc  G(L) 
u  o  o  o 


r  -°oco  fi(L>  i 

[2  +  p0c0  f i,(L)J 


where  G(L)0,  is  the  input  admittance  of  a  nonmagnetic  layer.  An  explicit 
expression  for  this  function  may  be  obtained  by  observing  that  according  to 
Equations  (46)  and  (52): 


fl(L) 


2i  j  sin  (kL  -  ®i)  sin(®2)  ) 
"  ”  pc  |  sin  (kL  “  "  «2)  J 


Substituting  in  Equation  (59)  and  using  Equations  (44)  and  (45),  we  obtain  the 
well-known  expressaion: 


(  G(0)  -  j 
ll  -  ipc  G 


yy  tan  (kL) 
G(0)  tan  (kL] 


With  the  aid  of  Equations  (58)  and  (59)  we  now  find  the  desired  expressions  for 
8,  80,  and  48: 


[1  -  P0c0  G(L)T 
1  ♦  ooco  C(L)J 


»  CU), 


where: 


r  1  “  Poco  C(L)o1 

1  ♦  P.c.  G(L) 
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n 


from  Equation  (59)  and: 


*i  r  t2(p  i 

-  *2  L2  *  Vo 


Hence  from  Equations  (59),  (61)  and  (63): 


Combining  Equations  (59)  and  (63)  we  also  obtain: 


48  f 
*o  “  y 


Poco  9l  i2^L^ 

»  c  £.(L))  (1 
o  O  1 


“  *2)} 


At  this  point  we  have  an  explicit  expression  for  R,  the  change  of  the 
reflection  coefficient  due  to  the  magnetization  Mq,  in  terms  of  the  Fourier 
transform  of  the  vector  potential  in  Equation  14.  Hence,  as  a  final  step  we 
calculate  the  quantity  Aa(z,z0,M)/a(e)  ■  ia(z,z0,«)/»Zo  from  the  gradient  of 
the  applied  potential  in  Equation  (11),  i.e.,  from  Aa0(z,z0,«)  “  a(u)  • 
8ao(s,so,«)/Szoiu*0.  Because  of  our  approxisiation  of  u(tQ,u)  it  is  not 
necessary  to  perform  the  differentiation  on  z0,  and  we  calculate  the  resultant 
potential  a(z,zQ,e)  from  the  applied  potential  a0(z,z0)e~i"t.  Clearly,  if  we 
did  not  want  to  use  the  average  of  u  in  the  dynamic  equations  we  could  obtain  the 
integrand  in  the  expression  for  the  magnetic  pressure  by  simply  differentiating 
the  expression  for  a(z,zOIu),  as  calculated  below,  with  respect  to  zQ,  and 
multiplying  by  u(z0,e).  Conceptually  the  potential  a(z,z0,u)  is  the  result  of 
applying  a  bifilar  current  at  z"z0  with  unit  amplitude  and  frequency  u. 

Physically  we  are  dealing  with  the  motion  of  a  constant  current  at  z  ■  z0  + 
u(z0,t).  This  current  is  changing  its  flux  linkage  with  the  ground  plane  and 
hence  is  inducing  eddy  currents.  If  the  motion  >u(z0,t)/8t  were  uniform 
we  could  have  derived  the  form  of  6a0(z,z0,t)  from  the  classical  theory  of 
electroawgnetism  in  moving  media.  For  if  the  applied  flux  is  denoted  by  6  we 
have  for  any  arbitrary  current  filament  in  the  plane  of  the  boundary  with  area 
element  da  and  line  element  ds: 


^  “  ft  fUo  *  d*  "  -  /E  o'  d® 

-  -  /?  X(v  X  b0)  •  n  da  -  -  / (v  X  bD)  •  ds  , 
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.  X  b0(z,z0) 


Aa0  “  ”  u  X  b0  *  -u  X  (boz  k  +  boy  j)  *  iu  boy 


since  a0,  is  a  function  of  (z-z0)  and  the  subscript  zero  on  the  electric 
field  Eq,  the  magnetic  field  b0  refers  to  "applied"  quantities.  Equation  (68) 
was  obtained  previously.  Equation  (11),  in  spite  of  the  fact  that  in  the  present 
case  motion  was  assumed  to  be  uniform.  This  agreement  is  evidently  due  to  the 
approximation  we  made  in  the  derivation  of  Equation  (11).  Had  we  gone  to  second 
order,  we  would  have  picked  up  a  term  of  the  form  lu(zQ,t)l2  92a0/9z02 iu“0 
which  would  have  introduced  the  curvature  of  the  potential  as  well  as  the  first 
harmonic,  with  frequency  2w. 

The  differential  equations  for  a(z,z0,w)  are: 


*  >  0}  F2a  »  -p0I04(z  -  zQ)  lA(y  -  A)  -  «(y  +  A)J 


z  £  0;  V2a  “  u0o  ~  ■  -iucrv0a  (70] 

where  the  time  dependence  e“lwt  has  been  cancelled  on  both  sides.  The 
Fourier  transform  of  the  source  potential  is  logarithmic  and  may  be  represented 
by  the  following  integral: 


l*o*o  C 

ao  J 


-k<z-z0l  sjn(kA)  sin(ky)  dk 


Solving  the  boundary  value  problem  by  requiring  the  continuity  of  a(z,z0,w) 
and  the  parallel  component  of  the  magnetic  field,  i.e.,  of  the  gradient 
9a(z,z0,«)/9z  at  z  *  o  leads  to  the  Fourier  transform  of  the  total 
potential  a(z,z0,m).  When  z  >,  o: 
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/  \  t  \  **0^0  C  'k(z+zo> 

a(z,z0,w)  ■  a(z,z0)  ♦  — —  /  e 


sinCkA)  ain(k; 
k 


1  (Hf)  -  ■ 


where  q  »  (k^  -  iwov0)*^.  In  passing,  we  note  that  the  corresponding 
solution  for  the  case  of  a  circular  disk  of  radius  is: 


y  ,  ^o^o^l 

a(z,z0,w)  - - 2 - 


/  {^'  * (Hf) 


-k(z+z. 


*  J^k^)  JjCkr)  dk 


where  J}(x)  is  the  cylindrical  Bessel  function  of  order  one  which  is  regular 
at  x  -  0. 

Note  that  the  absolute  magnitudes  of  the  integrands  decrease  exponentially 
with  increasing  k.  The  important  contributions  to  the  values  of  the  integrals 
therefore  occur  in  the  neighborhood  of  k  *  0.  The  intergals  are  readily 
evaluated  if  in  the  neighborhood  of  k  “  0  the  expression  ((k-q)/(k+q))  is 
approximated  by  an  exponential  function  with  an  exponent  linear  in  k.  We 
therefore  write: 


f(k)  •  ^  .  exp{g(k)] 


and  expand  g(k)  *  log  [f(k>]  in  a  Taylor  series  near  k“0  to  first  order  in  k. 
We  obtain: 


g(0)  «  logtf(O)]  -  0 


f 


f (k)  dk 


iwoyolk“0 


yS't  «*¥<:< 


Hence: 


f(k)  «  exp[-k(l  +  i)A]  ,  42  -  —  (75) 

UOUo 

and  6  is  the  electromagnetic  "skin  depth."  Substituting  this  into  a(z,z0,u) 
and  integrating  one  finds: 


a(z,z0,w) 


Mo  (!"(*  -  M2  +  (y  +  a)2 
*  08  (  L<z  -  z0)2  +  (y  -  A)2_ 


* 


(z  +  zQ  +  x)2  +  (y  +  A)2  ) 

(z  ♦  z0  +  x)2  ♦  (y  -  A)2j I  x  .  (1  +  iU 


(76) 


Evidently,  the  part  of  a(z,z0,o>)  in  Equation  (76)  that  is  independent  of 
4  is  due  to  internal  sources.  The  only  reason  for  carrying  it  to  this  point 
in  the  calculation  was  that  it  was  needed  as  the  source-term  for  the  eddy 
current  problem.  Its  effect  on  the  dynamic  response  of  the  plate  has  already 
been  accounted  for  in  Hooke's  Law,  Equation  (16a).  We  will  therefore  discard  it 
now. 


Substituting  the  remaining  term  in  Equation  (76)  into  Equation  (18)  and 
setting  y  *  +  A  (y  “  -A  was  accounted  for  when  the  volume  force  was 
calculated)  one  finds: 


P m(z,«) 


2I02au0  ( _ z  +  L  ♦  x _ 1 

A*  |(z  +  L  +  x)2  +  (2A)2  '  z  +  L  +  x 

(z  ♦  x)  +  1  ) 

(z  +  x)2  +  ( 2 A ) 2  z  +  xj  x  »  (i  +  i)s 


(77) 


where  |i0I02-B2*(gauss)/4v  is  the  magnetic  field  energy  density  per  cm^ 
in  cgs  units. 

The  expression  for  Pm(z,w)  may  now  be  substituted  into  Equation  (61)  to 
yield  the  reflection  coefficient  for  a  normally  incident  pressure  wave  and  into 
R,  Equation  (57),  to  yield  the  change  in  performance. 

We  have  evaluated  R  for  several  different  combinations  of  p,  c,  B  and 
A,  see  Figures  2  through  9.  The  plate  was  assumed  to  be  10  cm  thick  and  20  cm 
wide,  i.e.,  both  L  and  A  in  Figure  1  were  10  cm.  The  XFIL  integration 
routine  of  the  NSWC/WOL  computer  library  was  used  to  calculate  the  various 
integrals  in  Equation  (64). 


FIGURE  9.  ACOUSTIC  REFLECTION  FROM  VISCOELASTIC  AND  MAGNETOVISCOELASTIC 


A*  indicated  in  the  text  we  expected  to  have  to  use  rather  large  internal 
fields  in  order  to  obtain  observable  effects.  The  results  bear  this  out;  unless 
we  apply  at  least  several  hundred  thousand  gauss,  unless  we  use  an  Mq  that 
produces  a  B  of  this  order  of  magnitude,  there  are  no  discernible  differences  in 
the  responses  of  magnetic  and  nonmagnetic  plates.  We  find  that  the  energy  ratio 
in  Equation  (36)  is  indeed  the  controlling  parameter  in  this  situation.  The 
conductivity  of  the  substrate  is  of  relatively  minor  importance. 

In  closing  we  might  add  that  although  we  believe  the  order  of  magnitude  of 
the  effect  was  calculated  correctly,  we  could  have  done  better  had  we  used  the 
perturbation  expansion  of  the  solution  in  our  numerical  calculations.  In 
retrospect  we  feel  that  the  phase  information  that  was  suppressed  when  the  local 
displacement  in  the  magnetic  volume  force  was  averaged,  might  have  had  a 
significant  influence  on  the  details  of  the  frequency  dependence  of  B^.  On 
the  other  hand,  in  view  of  the  many  conceptual  simplifications  we  have  made  in 
the  course  of  solving  this  problem,  this  refinement  would  have  been  academic. 
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CHAPTER  2 

ACOUSTIC  ENERGY  ABSORPTION  IN  AN  ELECTRICALLY  CONDUCTING,  VISCOELASTIC 
MATERIAL  IN  A  UNIFORM  MAGNETIC  FIELD 


In  Chapter  I  the  sources  of  the  tine  varying  magnetic  field  on  the  one  hand 
and  the  electromagnetic  absorber  on  the  other  were  spatially  separate.  If  the 
aagnetoviscoelastic  material  was  made  electrically  conducting  both  functions 
could  be  combined  in  the  same  volume.  Two  apparently  different  mechanisms  would 
now  give  rise  to  conduction  currents.  In  addition  to  the  mechanism  discussed  in 
Chapter  1,  which  we  think  of  as  a  9B/3t  effect,  there  would  now  be  the  (v  X 
B)  effect,  which  is  due  to  convective  transport  of  charge.  Both  effects  are  due 
to  the  relative  motion  of  "sources"  and  "observers."  In  the  first  case  the 
"sources"  awe,  hence  SB/ at  at  the  observer  was  picked  up  as  (v  X  B)  in 
Chapter  1,  where  v  referred  to  the  motion  of  the  source  (the  magnetic  moment). 

In  the  second  case  the  "observer"  moves,  i.e.,  v  refers  to  the  motion  of  the 
observer  (the  mobile  charge  carrier).  Note  that  in  the  former  case  relative 
motion  of  source  and  observer  is  detectable  only  because  B  is  nonuniform. 

The  method  of  analysis  used  in  Chapter  1  is  too  approximate  for  a  situation 
in  which  the  awgnetic  material  discussed  in  that  section  simply  becomes 
conducting.  The  direction  of  the  magnetic  field  inside  the  magnetized  slab 
would  also  be  wrong,  since  it  is  roughly  in  the  same  direction  as  the  local 
displacement  and  therefore  eliminates  the  (v  X  B)  effect..  If  on  the  other  hand 
the  uniform  magnetization  is  assumed  to  be  parallel  to  the  ground  plane  instead 
of  perpendicular  to  it,  then  the  internal  induction  will  also  be  roughly  uniform 
and  parallel  to  the  ground  plane  and  effects  due  to  finite  conductivity  are 
maximised.  These  effects  may  then  be  analyzed  readily  if  at  the  same  time  we 
neglect  the  magnetic  field  variations  that  were  the  source  of  dissipation  in 
Chapter  1,  namely  those  due  to  local  displacements  of  magnetic  moments  under  the 
influence  of  dilatational  waves  and  focus  only  on  convective  charge  transport 
arising  from  the  translational  motion  of  the  molecules;  i.e.,  on  the  (v  X  B) 
effect.  The  constant  swgnetization  M0  may  then  be  replaced  by  the  constant 
amgnetic  induction,  BQ,  where  B0  *  m0Mq,  or  equivalently  we  may 
consider  an  unmagnetized  material  in  a  constant  parallel  field  with  a  wave 
travelling  at  right  angles  to  the  field  direction.  In  fact,  since  we  really 
only  need  the  attenuation  constant,  or  the  imaginary  part  of  the  propagation 
vector,  k,  in  order  to  estimate  changes  in  power  dissipation,  we  will  only  solve 
for  a  wave  in  an  infinite  medium  rather  than  for  reflections  from  a  layer. 

Hence* we  are  looking  for  the  complex  dispersion  relation  k(u,o,B). 

Under  the  action  of  a  pressure  wave  in  the  material  we  think  of  the  mobile 
charge  as  being  "carried  along"  by  local  displacements  of  molecules,  i.e.,  there 
is  no  charge  separation.  The  (v  X  B)  electric  field  creates  a  current,  J,  and 
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this  current  locally  exerts  a  Lorentz  force,  J  X  B  on  the  material.  As  in 
Chapter  1  the  equations  of  motion  are*. 


.£--£*»uxhi 


(78) 


and: 


fiC 


2  H  .  .  !£ 


iz 


at 


(79) 


As  in  the  theory  of  magnetohydrodynamics,  the  equation  for  the  magnetic 
field  is  obtained  from: 


T  X  1  *  |iJ 


(80) 

(81) 


and: 


J  -  o[E  ♦  (v  X  B)] 


(82) 


in  aka  units,  where  the  last  equation  is  stricly  true  only  for  a  uniform 
velocity,  v.  Eliminating  J  and  E  from  Equations  (80),  (81),  and  (82)  yields: 


~  ?2B  ♦  V  X  (v  X  B)  -  ||  (83 

for  the  magnetic  induction.  He  assume  that  B  consists  of  a  strong,  stationary 
background  field,  B?,  and  a  relatively  weak  response  field,  b;  B  *  B0  +  b. 
Furthermore,  if  BQ  is  assumed  to  point  in  the  x-direction  and  v  in  the 
z~direction,  then  J  is  in  the  y-direction  and  the  response  field,  b  is  parallel 
to  Bc. 

Assuming. the  time  dependence  e“*-at  in  Equations  (78),  (79),  and  (83) 
yields: 


iw(F  X  H)  X  B  -  pc2 


dv 

dz 


pw2v 


(84  j 
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and: 


~  72B  ♦  V  X  (v  X  B)  -  -iwb 
t>a 


(85 


Taking  account  o£  the  fact  that  the  only  coordinate  dependence  is  on  z,  we  find: 


-in  d(B2)  ^  d2v  +  id2 
2ppc2  dz  dz2  c2 


and: 


(86; 


po 


d(vB) 

dz 


— r  +  iwpab 
dz2 


(87 ; 


Using  the  fact  that  BQ,  is  constant,  Ib/B^l  «  1  and  neglecting  products 
involving  v  and  b  compared  with  products  involving  BQ,  in  combination  with  v 
or  b,  we  obtain: 


d2b 

dz2 


♦ 


iuyob 


(88) 


(89) 


The  desired  dispersion  relation  for  the  propagation  of  dilatational  waves 
in  the  material  is  now  obtained  if  we  substitute  the  solution. 


v  *  vQ  e**2  ,  b  *  bQ  e*-kz  (90 

into  Equations  (88)  and  (89): 

[*  -  g]  p  -  H  -  ^  - 

or: 

lw(  o  v  V2 

k  •  ♦  ” TF’T  U1  +  i  x(l  ♦  y))  +  [(1  +  i  x(l  +  y))Z  -  Aix]1/2}  (92 
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riiere: 


HoBo 

y  -  — Y~ 
pc * 


The  expression  in  the  curly  bracket  of  Equation  (92)  therefore  modifies  the 
propagation  constant  k0  *  w/c  which  holds  in  the  absence  of  a  background 
field  or  when  0*0.  Note  that  the  parameter  x  may  also  be  written: 


where: 


5  ■  V2 /uo» 


is  the  electromagnetic  "skin  depth,"  or: 


-  (tf 


x  ■  5  x  10 


where  XQ  is  the  wavelength  in  the  material  in  the  absence  of  magnetic  effects. 
The  two  parameters,  y  and  4  are  seen  to  be  the  controlling  variables  of  this 
problem,  as  they  were  in  the  reflectivity  problem  of  section  (A).  Let: 


4v  Henri 


a  ■  102 

m 


n  *  2i  x  10*  Hz 


1.5  «  103 


Bc  -  10*  gauss  «  10'1 


p  -  i 


10* 


Then: 


y  -  3.54  x  10“6 

Hence. with  these  orders  of  magnitudes  for  the  parameters,  the  ratio  of  magneto¬ 
static  to  kinetic  energy  density  is  very  small;  y**o  and  from  Equation  (92)  it 
folows  that  the  dispersion  relation  is  the  same  as  in  the  case  of  or  *  o. 

k  ~  H- 
c 

The  induction  inside  the  material  has  to  rise  to  the  order  of  magnitude  of 
several  hundred  thousand  gauss  before  the  dispersion  relation,  and  hence  the 
absorption  mechanisms  are  significantly  altered.  p0te  also  from  Equation  (88) 
that : 


Since  vQ  in  this  type  of  situation  is  of  the  order  of  magnitude  of  10'  to 
10“^  m/sec,  the  assumption  that  l b/B0 I  «  1  is  satisfied.  Evidently,  the 
inequality  remains  valid  even  when  x  <<  1. 

In  Figures  10  through  18  we  have  plotted  the  real  and  imaginary  parts  of  k, 
Equation  (92),  over  the  real  and  imaginary  parts  of  (u/c)  respectively  as 
functions  of  frequency. 

As  in  Chapter  1  we  find  that  large  magnetic  fields  are  needed  for  strong 
effects. 


Summarising  the  results  of  both  chapters  we  conclude  that  unless  the  ratio 
of  the  characteristic  magnetic  field  energy  density  associated  with  a  given 
situation  to  the  mechanical  energy  density,  pc^,  of  the  material  is  at  least 
of  order  unity,  neither  the  electrical  conductivity  of  a  magnetoviscoelastic 
material  nor  induced  eddy  currents  in  neighboring  structures  should  be  expected 
to  lead  to  significant  electromagnetic  absorption  of  elastic  waves  or 
dissipation  of  acoustic  energy. 


FIGURE  10.  OISPERSION  RELATION  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  (Bft“  105  GAUSS,  a  -  104  MHO/M,  pc2/pnc£  -  1AMl.-W.1i.  p/p0  -  5.C 


2  (0)  coxjui/ooui  (*)  (oxnyaa/oi) iwaa 
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FIGURE  11.  DISPERSION  RELATION  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  (B0-  105  GAUSS,  o  -  104  MHO/M,  pcZ/p0cg -1.0* (1 .10.1).  p/p0  -  1.01 


FIGURE  12.  DISFERSION  RELATION  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  (Bg“  106  GAUSS,  o  -  104  MHO/M,  pe2/pflc§  -  0.1  Ml. -MM).  p/p0  -  1.0) 


FIGURE  14.  DISPERSION  RELATIONS  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 

CONDUCTOR  (B0«  S  x  105  GAUSS,  a  -  104  MHO/M,  pe2/p0e^  -  0.1*(1.-»0.1>.  plpQ  -  1.0) 


FIGURE  15.  DISPERSION  RELATIONS  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  (B0-  10®  GAUSS,  o  -  105  MHO/M,  pc2/p0cg  -  0.1  *(1.-10.11.  p/pQ  -  1.01 


FIGURE  16.  DISPERSION  RELATIONS  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  (Bn«  5  x  105  GAUSS,  o  -  105  MHO/M,  pc2/p„eS  -  0.1*11.-10.1 ).  plpn  - 
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FIGURE  17.  DISPERSION  RELATIONS  FOR  SOUND  WAVES  IN  A  VISCOELASTIC  ELECTRICAL 
CONDUCTOR  |B0-  10®  GAUSS,  o  -  10®  MHO/M,  pc2/p„cg  -  0.2*11.-10.1 ),  p/pQ  -  1.0) 
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